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We discuss in depth the validity and limitations of a theoretical scheme to evaluate photo-electron 
spectra (PES) through collecting the phase oscillations at a given measuring point. Problems appear 
if the laser pulse is still active when the first bunches of outgoing flow reach the measuring point. This 
limits the simple scheme for evaluation of PES to low and moderate laser intensities. Using a model 
of free particle plus dipole field, we develop a generalized scheme which is shown to considerably 
improve the results for high intensities. 

PACS numbers: 33.60.+q, 33. 80. Eh, 36.20.Kd 
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I. INTRODUCTION 



Photo-electron spectra (PES) constitute a longstand- 
ing basic tool for analyzing the electronic structure of 
atoms, molecules, or solids [JQ- In the one-photon 
regime, PES provide an image of the sequence of single- 
particle levels which are occupied in the ground state. 
This has been applied, e.g., already in the early days of 
cluster physics to the electronic structure of cluster an- 
ions to track the transition to bulk metal [||. While this 
could be done with photon frequencies in the range of 
visible light, the analysis of deeper levels requires higher 
frequencies, and hence one finds also UV [5[ and X-ray 
PES [f|. Nowadays, we dispose of a great variety of co- 
herent light sources in large ranges of frequency, intensity 
and pulse length, such as the very powerful and versa- 
tile free electron lasers @, HI which now even allow for 
time-resolved studies of deep lying core states of atoms 
embedded in a material @. With the great availability 
of good light sources, studies of PES are now found in 
all areas of molecular physics, from atoms over simple 
molecules [l(| to complex systems as clusters [llj and 
organic molecules [T3 |. 

There is hence a general need for a robust theoretical tool 
of analysis of PES in various dynamical regimes. Tradi- 
tional approaches to compute PES rely on (multi-)photon 
perturbation theory (l3| . A few years ago, we have de- 
veloped a technique to evaluate the PES directly from 
numerical simulations of the electronic excitations on a 
spatial grid representation (and also for a plane- wave rep- 
resentation) [lUdH allowing exploration of a rather wide 
range of dynamical scenarios [lj] . The basic idea of the 
procedure is to use absorbing boundary conditions and 
to assume free particle propagation near the boundaries. 
The technique has indeed been successfully applied to a 
variety of clusters 0, [TtJ , for which significant electron 
emission can be achieved with comparatively low laser 
intensities. Atoms and molecules with larger ionization 
potentials (IP), as e.g. carbon, require stronger fields in- 
stead and it was found in first tests that the application of 
the recipe from Ref. [l4| can yield artefacts for the case of 



long pulses with high intensity. Therefore it is a timely 
task to re-inspect the method and to develop improve- 
ments where found necessary. Such a critical analysis of 
the evaluation of PES on a numerical grid is the aim of 
this paper. 

We start in section [IT] with a brief review of the tra- 
ditional recipe to evaluate PES. We show that, while 
this simple recipe works well for low laser intensities, it 
produces unphysical spectral features at high intensities, 
which call for a revision of the method. We continue 
with a quick reminder of the gauge freedom in describ- 
ing the laser field. We discuss two choices of gauge and 
explain the gauge transformation that connects the two 
choices. The question which gauge is ideally to be used 
in the evaluation of PES will be answered in the course 
of the further sections. Afterwards, the presence of a 
strong laser field in coexistence with the emitted electrons 
is dealt with by solving the time-dependent Schrodinger 
equation for free electrons plus a homogeneous laser field. 
This case still allows closed solutions which is then used 
as a basic ingredient to define a refined recipe to evaluate 
PES. In section IITT1 the refined recipe is then tested for 
the analytically solvable case of Gaussian wave packets 
first and then numerically for a few realistic test cases. 



II. THEORETICAL EVALUATION OF PES 



In this section, we first give a brief review of the scheme 
for the evaluation of PES introduced in [3, EH ■ After 
showing that it is not appropriate for high laser inten- 
sities, we propose a generalized scheme that solves the 
problem. All schemes discussed here and in the previ- 
ous papers rely on a mean-field description where each 
electron is associated with a single-electron wave func- 
tion ip a . We will consider in the formal discussions one 
representative state and drop the index a to simplify no- 
tations. 
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A. PES scheme for free electron 

The PES are calculated in a very efficient manner ex- 
ploiting the features of the absorbing boundary condi- 
tions [Til . [lH . We give here a brief summary of the pro- 
cedure and its motivation, first for the simple case of ID 
and then in 3D. It is to be noted that these traditional 
recipes for evaluating the PES arc deduced under the as- 
sumption of free particle propagation near the bounds of 
the numerical box. 



1. The ID case 

We choose a "measuring point" zm far away from the 
center and one or two grid points before the absorbing 
boundaries. We record the wave function ip(z,t) at the 
measuring point zm during the time evolution. This 
delivers the raw information from which we extract the 
PES. In order to develop an appropriate recipe, we now 
need a few formal considerations. 

We assume that the mean-field is negligible at zm such 
that we encounter there free particle dynamics which is 
governed by the one-particle time-dependent Schrodinger 
equation idtip(z,t) = —^-^p(z 7 t). The solution is 



i/>(z,t) 



dk 



(kz—ujt) 



(1) 



with the condition that k and uj satisfy the dispersion re- 
lation u> = k 2 /2. Note that atomic units (h = m e = e = 
4-7T60 = 1) are used here and throughout the paper. The 
measurement of PES is practically a momentum analysis 
of the outgoing wave packet at a remote side. In momen- 
tum space, the wave function (fTJ) reads 



'2-k 



e~ lH:x ip(z,t) = tP(k,0)e~ 

ocg(k) 



(2) 



The probability to find an outgoing particle with momen- 
tum k is thus given by \i/j{k,t)\ 2 = |V>(fc,0)| 2 , from which 
the PES can be obtained as 



y(E kin ) cx 



V(fc,o) 



(3) 



where -Ekin = k 2 /2 = u> is the kinetic energy and we 
have taken into account the appropriate energy density 

— 1/2 

cx E kir ^ . This is what we call in this work the exact 
definition of PES developed close to the experimental 
procedure. 

However, it is numerically extremely expensive to eval- 
uate the PES by Fourier transformation for a remote slot 
in coordinate space. To overcome this difficulty, we ob- 
tain the expansion coefficient g(k) from the Fourier trans- 
form in time- frequency space of the wave function ij)(z, i) 
collected at a measuring point zm near the absorbing 



boundary. We denote the Fourier transformation in time 
by -0 to distinguish it from the spatial Fourier transform 
i/>. This then reads : 



(1/ 



2tt 
At 



ip(zM,t) e 



iui t 



/2tt 
dk 



dk 
/2n 



9{k) 



,i(kzM—ut) iu>' t 



g(k) e ikzM 5{oj-u)') 



(4) 



Since we collect ip(zM,t) far away from the center and 
close to the absorbing boundary, we can assume that only 
outgoing waves will pass the measuring point zm- We 
will hence have only positive wave vectors 
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(5) 



which allows us to approximate the last line in (j4]) as 

du 



ip(z M ,u') 



■.g(uj)e lZMV2LJ S(lj - uf) 



TTUJ' 



(6) 



In the last line, 9 stands for the Heaviside function. Using 
© and (@|, we finally obtain 



-i zm \/2uj 



ip(z,uj) 



(7) 



where from now on, we will always consider u > and 
thus omit the factor 0(u>) in Eq. ([7]). At measuring point 

zm, we can identify ip(zM,u) with ip{zMi -E-kin) where 
i?kin is the kinetic energy of the electron. This defines 
the PES yield y ZM (E kin ) as 



yz M (E kin ) OC V^kin -ipizMiEkin) 



(8) 



We have checked the method in extensive ID wave packet 
calculations and compared it to a direct momentum de- 
composition of the outgoing wave as given in Eq. ([3]). 
Both methods yield the same results, while the above 
sketched frequency analysis at a " measuring point" is or- 
ders of magnitude faster. 



2. The 3D case 

As in ID, we assume for the 3D case free propaga- 
tion and purely outgoing waves at the measuring point 
ym- The solution of the one-particle timc-dcpcndcnt 
Schrodinger equation in free space reads 



ff>(T,t) 



d 3 k 



ff(k)e 



i(k-r— tot) 



(9) 



Analogously as in Eq. ([5]) for the ID case, we calculate 
the coefficients g(k) from the time Fourier transform of 
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ip(r,t), assuming that only wave vectors k = ke r with 
k > contribute, where e r is the direction of the outgoing 
radial wave. This yields 



(10) 



where f2 r is the solid angle related to the r direction. 
Considering ym as the measuring point, we compute the 
PES in direction £l r .. as 



yn IM (Skin) oc ip(r M ,E kin ) 



(11) 



The above analysis yields the fully energy- and 
angular- resolved PES. The angular averaged PES would 
then be attained by angular integration with proper solid 
angle weights. It has nevertheless to be emphasized 
that the resulting angular dependence in the laboratory 
frame depends on the orientation of the molecule (or clus- 
ter). Actual ensembles in the gas phase do not contain 
molecules in well defined orientation but represent rather 
an equi-distribution of all possible orientations. A typical 
example are the many recent measurements of angular- 
resolved PES in cluster physics see, e.g., fl8l - [22l ]. An 
appropriate orientation averaging has to be performed 
before one can compare yn IM (Skin) with experimental 
data. The practical procedures for that are outlined in 
[23l l24j . Orientation averaging, however, is beyond the 
scope of the present paper and will be ignored in the 
following. 
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3. Example of application and problem 

Although this simple scheme has been applied with 
success to a variety of clusters, our recent attempts to 
compute the PES of small covalent molecules have raised 
some questions concerning its general applicability. For 
high laser intensities, we find an unexpected shoulder in 
the PES at large kinetic energies. These pattern were 
then also found for the simple test case of Na n clus- 
ters with jellium background when going to sufficiently 
large laser field strengths. This excludes particularities of 
pseudo-potentials, local or non-local ones. The problem 
seems to reside in the scheme to evaluate PES. 

FigureQ]dcmonstratcs the problem for the case of Na^~. 
The spherical jellium model is used for the ionic back- 
ground. Valence electrons are described by the time- 
dependent local-density approximation (TDLDA) using 
the energy functional of |25[ and an average-density self- 
interaction correction (ADSIC) (2(|. Two cases are con- 
sidered for comparison: one with still moderate laser in- 
tensity, and another one in the high intensity regime. In 
the latter regime, more than half of the cluster's electrons 
are stripped off, as shown in the bottom panel of Figure 
[TJ The top panel shows the PES for both cases. For mod- 
erate intensity, we see the typical monotonous decrease. 
The case of high intensity differs : Up to Sdn « 1.5 Ry, 
we see the typical pattern of a monotonous decrease of 



FIG. 1: Ionization properties of Nag with jellium back- 
ground under the influence of laser pulses having frequency 
wias = 0.1 Ry, pulse length T pu i se = 12 fs, and intensities 
7 = 10 13 W/cm 2 (green or light gray lines) and I = 10 14 
W/cm 2 (brown or dark gray lines), computed in a cylindrical 
box of 176x88 a 2 , with spherical absorbing bounds covering 
at least 16 grid points. The laser field was given effectively 
in u-gauge. Bottom : Time evolution of ionization. Left : 
Photoelectron spectra y» TM (Skin), see Eq. |TT|). 



the envelope with some fluctuations. For E^ n > 1-8 Ry, 
however, we observe a new maximum, a broad shoulder 
of high energy electrons; this is totally unexpected and 
most probably unphysical. In the following, we clarify the 
origin of such a shoulder and we generalize the method 
to compute PES to a wider range of laser intensities. 



B. Space gauge and velocity gauge 

Before proceeding with a cure for the PES scheme illus- 
trated above, we review the choices of gauge for describ- 
ing the laser field. The external laser field is described 
in the limit of long wave lengths and neglecting mag- 
netic effects (electronic velocities are very small at any 
time). The electron- laser interaction can be described 
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either by a time-dependent scalar potential $(r, t) or a 
time-dependent vector potential A(r, t) in the Hamilto- 
nian, which in general reads 



H(v,t) = 



1 



A(r,t) 



i 2 



(12) 



In the space gauge (x-gaugc) , the laser field is described, 
within the dipole approximation, by a scalar potential 
only: 



A<*>(r,t)=0, &*\T,t) = -Eof{t)i 



(13a) 



where the laser polarization is chosen along the z axis, 
and Eq is the maximum electric field of the laser. f(t) 
here denotes the temporal profile of the laser pulse, taken 



f(t) = sin(w las £) sin 2 



7T t 
^pulsc 



9(t) 6(T ulse -t) , (13b) 



where wi as and T pu i so are respectively the frequency and 
the duration of the laser pulse. 9 denotes the Heaviside 
function. We have chosen a smooth sin 2 profile which 
combines high spectral selectivity and finite extension. 
The spectral selectivity is needed because the pulse pro- 
file carries through to the PES 15 1. 

In the velocity gauge (w-gauge) instead, the laser field 
enters the Hamiltonian as a vector potential only: 



A("> (r, t) = -cE F(t) e z , (r, t) = , (14a) 



with 



F(t) 



dt'f(t>). 



(14b) 



Note that the vector potential in dipole approximation is 
spatially constant. 

Both gauges are equivalent. They are connected by a 
gauge transformation which in general reads : 



X 



A = A-V X , 



ip = ip exp 



(15a) 
(15b) 
(15c) 



Let us assume that we take the w-gauge as a starting 
point : the laser field is described by the vector poten- 
tial AM = -cE F(t)e x , whereas $W = 0. If we want 
to gauge transform the vector potential into the scalar 
potential, then Vx = A" , from where 

X = -cE F(t)z => <f>^ = ^ = -E f(t)z . (16) 

c 

This finally gives 

</,<» = ,/,(«) exp [- iE F(t)z] (17) 



which relates the wave function tp^ as computed by 
time propagation under the action of the scalar poten- 
tial (|13a|) . and ip( v > for the case with the vector potential 
(|T4"|). It is obvious that this phase factor is crucial in 
the present PES analysis, while it can easily be ignored 
for local observables as, e.g., dipole momentum or net 
ionization. 

The phase transformation allows one to decouple the 
gauge used in the analysis from that used in the time 
evolution. Assume that we solve the Schrodinger equa- 
tion using the potentials (|14|) . This immediately yields 
the wave function ip^ v > in u-gauge. The same wave func- 
tion could be obtained by using (fT3|) and applying the 
reverse transformation (|T7|) to recover om the tpW 

as obtained by propagation. In fact, this is the most effi- 
cient way to evaluate tp^ as the operator (|13ap is purely 
local. It is to be noted that the transformation (fl7|) is 
relevant for the PES only if the outgoing wave reaches 
the measuring point zm at a time where the signal F(t) 
is still active. For very large boxes, the outgoing wave 
and the signal can avoid to coincide, in which case each 
gauge yields the same result. However, typical box sizes 
used in practice are not always that large. 

The question is now which gauge is most suitable for 
evaluating PES. The basic papers [13, H3 state that the 
preferable choice is the w-gauge. The argument there was 
a better convergence with box size. We now speculate 
that this has to do with (non-) coincidence of outgoing 
wave packet and laser signal. In the next section, this 
will be scrutinized using a solvable model. 



C. PES scheme for free particles plus laser field 

In section HH we have deduced the recipe for evaluating 
PES under the assumption that the potential is negligible 
at the measuring point. Although this may hold for the 
typical mean field of a system, we cannot easily exclude 
the presence of the laser field at this point because the 
laser field is of extremely long range (wave length much 
larger than system size in the dipole approximation that 
we use here). Thus we extend the considerations to the 
case of a free particle plus laser field in dipole approx- 
imation which is still analytically solvable. We confine 
the considerations to the ID case. The extension to 3D 
is straightforward. 



1. Momentum-space wave function in v-gauge 

We first consider the Schrodinger equation of a free 
particle with the external field ([H]) in v-gauge 



l[p-E F(t)e z ] 2 ^(z,t) = id t ^(z,t) ,(18a) 



(18b) 



5 



It is advantageous to expand the wave function at a given 
time in momentum space: 



V> (u) (M) = / -^Le~ ik * ^\z,t) 



with 



^ (t,) (fc,0) = fy v \k) . (19b) 
This yields the Schrodinger equation in the form 

X - [ft - E F(t)} 2 (ft, f ) = i (ft, t) , (20) 

which solution reads : 

(ft, t) = cxp (- i jf dt' [ft - EoF(i')] 2 ) ffi (ft) • 

(21) 

We rewrite the latter expression as 

^ [v) (ft, t) = cxp ( - iut + ikSq - iSfl) (ft) (22a) 
with 

" = y , (22b) 

^(i) = #o / dt'F(t') , (22c) 
Jo 

= ^ f dt'F{t') 2 . (22d) 
2 Jo 



Eqs. ([22]) constitute the solution when propagating the 
wave packet in v-gauge. 



2. Momentum- space wave function in x-gauge 

The Schrodinger equation of the same situation in x- 
gauge reads 



y + E zf(t) 



^ {x \z,t) = id t ^ x \z,t) ,(23a) 



yj {x} (z,0) = ^\z) , (23b) 



/ = dtF 



(23c) 



Note that at t = 0, the laser field is off and therefore, the 
wave function in x-gauge is the same as that in u-gaugc. 
In the following, we hence indicate the wave function in ft 
space at t = indifferently as V'o(ft) since the upper script 
(x) or (v) is irrelevant. Thanks to the gauge transforma- 
tion l|17p. the solution of (|23a[) can be easily obtained 
from the solution (|22ap in w-gauge, that is : 



ft x \k,t)=^(k + E F(t),t) . 
The wave function in x-gauge thus becomes 
$ {x) (k, t) = cxp (-iut + \K5q - Wfi) $ Q (k), 



(24) 



(25a) 



K = k + E a F(t) 
K 2 

U) = — 

2 ' 



(25b) 
(25c) 



where Sq and Sfl as given in Eqs. (|22c[) and (|22d[) re- 
spectively. This form is considerably more involved than 
the solution (|22a[) in u-gauge, since K explicitly depends 
on time t. We hence take this as a first formal indication 
that the preferred gauge to evaluate PES should be the 
w-gauge. 



3. Reconstruction of the PES 

The scheme to evaluate PES, as sketched in section 
ITU collects a wave function at a certain coordinate-space 
point ZjVj and performs Fourier transformation in fre- 
quency. This yields, by virtue of u; = fc 2 /2, the outgoing 
wave function in momentum space from which we then 
deduce the PES. In extension of the simple form ©, we 
now want to take care of the possible coincidence of the 
laser field at the measuring point with the by-passing 
outgoing wave. This setup has been solved in section 
III C II We are now going to apply the PES analysis to 
the solution (|22ap and deduce how to modify the collected 
signal at zm in order to properly describe the underlying 
unperturbed wave packet ipo(k). 

The relation between the wave function in frequency 
space and that in momentum space is given by 



dt 
dfc 



dt 



/2tt 



ty(M). (26) 



In the case where the laser field at the measuring point 
is negligible, one can see from Eqs. (|22a|) and (|22b>[) that 
the solution in fc-space shrinks to a simple ^ v '(k,t) = 
cxp (— ioj'i) ^o(ft) with u>' = Then the integration 
over t in (|26[) immediately produces a S(lu — a/), and that 
over ft, using the fact that ft = \J2uj' close to the absorb- 
ing boundary, yields ^ v '(z,u>) cx ipo(k). Therefore, the 
PES in Eq. © can be directly evaluated from ijj(z,u)) 
(see Eq. ©). 

Now we proceed to the case with non- vanishing laser 
field at the measuring point, once again in w-gauge. Try- 
ing to apply the time-frequency transformation directly 
to the fc-space solution (|22a|) runs into trouble due to 
the non-trivial time dependences induced by the factors 
Sq(t) and 8Cl(t). An obvious solution is simply to counter- 
weight the disturbing phase factors by a phase-correction 
factor exp (— ikSq + iSCl). This can be done even in coor- 
dinate space because these factors do not depend on the 
position. Therefore we calculate the PES from a phase- 
augmented Fourier transform ?/»( pA ), which is now pro- 



6 



portional to 



1. Gaussian wave packet in a laser field 



il>W(z M ,w) oc / dte iut e iv ^(z M ,t) , (27a) 
Lp = -kSq + SQ . (27b) 

This modified Fourier transform is then used as in Eqs. 
([8]) yielding the recipe for the ID case 



(28) 



The generalization for the 3D case proceeds analogously. 

The new move in this generalized evaluation is to aug- 
ment the original recipe by the phase-correction factor 
e lv . This should allow the application of the recipe in 
a wider range of laser intensities and time profiles. The 
factor becomes negligible for weak laser field, or for fields 
which do not interfere in time with the emitted parti- 
cle flow. We shall call in the following the former recipe 
© the "raw" recipe and the generalized form (|28p the 
"phase- augmented" (PA) recipe. A word of caution is 
in order. The additional phase involves the mere mo- 
mentum k. The recipe (|27a[) identifies that with the fre- 
quency u> according to Eq. ([5]). This is just as valid as 
the Fourier transformation is selective in frequency. This 
may be at stake in cases of very violent dynamics. 

We run into unsurmountablc trouble if we try to de- 
velop a generalized scheme on the basis of the solution 
()25a[) in x-gaugc because here even the instantaneous mo- 
mentum K, as given in Eq. (|25b[) . depends on time. This 
is another strong indication that the u-gaugc is the ap- 
propriate starting point for evaluating the PES (as was 
already argued in [l4|, H3] ) • We will confirm this suspicion 
in the next section with a practical example. 



III. ILLUSTRATION 

In this section, we show how the generalized scheme 
(|27a[) performs. We first apply it to an exactly solvable 
model, and then to more realistic cases. 



We consider a one dimensional system with a Hamilto- 
nian consisting of the free kinetic energy plus the external 
laser field, as it was discussed in section IPTCl The wave 
function of the system at initial time is taken as the fol- 
lowing Gaussian function: 



if ^ 1 (■ ( z -<io) 2 . n 

^ 0) = j^F 1 exp v poz " in ° 



(29) 



The Fourier transform in k space of (|29p reads 



Mk)=(^) 1 ' i e 



(30) 

The solution of the time-dependent Schrodingcr equation 
is detailed in Appendix [Al using the x-gaugc. It reads: 



Mo 



i/'-i 



exp 



[z ~ q(t)} 2 



o i(p -E F(t))z 



2(mo+«) 



with 



q(t) = q +Pot-E Q / dt'F(i') 



ifi(i) , (31) 



(32) 



n(t) = n Q + ^t- Po E Q J dt'F(t') 

+^! /*dt'F 2 (t') + iatanf— ) (33) 
2 Jo 2 \Mo/ 

Using now the gauge transformation (|17p , the wave func- 
tion in w-gauge becomes 



^\z,t) 



1/4 

Mo \ ip oZ 



7T(/i 2 + i 2 ) 



exp 



[z-q(t)\ 2 



2(mo + it) 



^- - ifi(t) . (34) 



A. Analytical test case: Gaussian wave packets 

As a first test case, we consider a Gaussian wave packet 
whose propagation can be analytically described, even in 
the presence of a (homogeneous) laser field. The ana- 
lytical solution, carried forth to far distance and beyond 
the lifetime of the laser field, allows a direct evaluation 
of the PES by filtering the momentum components by 
Fourier transformation from coordinate to momentum 
space. This "direct" analysis exactly corresponds to the 
experimental procedure, see Eq. ©. Thus we can test 
the PES analysis in frequency space against these "exact" 
results. 



2. Exact evaluation of PES by spatial Fourier analysis 

An exact evaluation of the PES has to analyze the 
momentum content of a wave packet at a place where 
free propagation is reached. The momentum distribution 
computed from the wave packet tp^ by spatial Fourier 
transformation becomes 

t) = (^) 1/4 e i-s(M)/2 exp (-|[fc-p(t)] 2 ) 
exp (i[p(t) - k]q(t) - ifi(t)) (35) 

where arg(/i) denotes the argument of [i. The momentum 
distribution and corresponding probability becomes in v- 
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gauge 



fi v) (k,t) = ( 



1/4 



e iarg(M)/2 eXP (-|[fc 



\PoJ 



f(J>o 



exp 

1/2 



(i[po-*!]?(t)-in(t)) , 



— ) exp (-fio[k -pa]' 



(36) 



The corresponding kinetic energy distribution is obtained 
by the identification k = +v / 2£m^. One also has to ac 

^ In 

count for the appropriate energy density oc i? kin 
yields 



- 1/2 This 



y [v \E Mn ) cx 



1 /2 

V0\ 1 e -MV2E^~-pof (37) 



This is then the exact energy distribution evaluated from 
spatial Fourier transform. 



- 



!Z3 

o. 



I=10 10 W/cm 2 




FIG. 2: Photoelectron spectra from wave packet analysis for a 
weak intensity I = 10 10 W/cm 2 . We use the "raw" recipe Q 
and compare analysis from the wave function in v- (brown or 
dark curve) and x-gauge (green or light gray line) . The exact 
result, that is the momentum distribution of the wave packet 
as given in Eq. (|37|) is also shown (black line). 



3. Evaluation with the phase-augmented PES scheme 

The final PES analysis for the above model relies on the 
wave functions at the measuring point, i.e. tp^ x '(zj^i,t) 
and i/j( v '(zM,t) as given in Eqs. and (pJTj) . respec- 

tively. We use the laser profile (|13b[) . As an actual ex- 
ample, we consider the following parameters: po = 1&q , 
go = 0, il = 0, uj las = 0.11 Ry = 1.4 eV, T pu i sc = 
8000 Ry" 1 = 384 fs, fi = 100 a and various field 
strengths Eq. The measuring point was taken at zm = 
1000 ao or 2000 a . The first choice yields an overlap be- 
tween laser pulse and flow signal at z^i while the second 
choice decouples them. In the following discussion, we 
quantify the field strength in terms of the laser intensity. 
The correspondence is, e.g., I = 10 14 W/cm 2 Eq = 
0.109 Ry/a , I = 10 12 W/cm 2 o E a = 0.0109 Ry/a , 
I = 10 10 W/cm 2 <-> E a = 0.00109 Ry/a - 

Figure [2] compares the evaluation of PES for x-gauge 
versus w-gauge for a case of weak laser field where the 
phase correction ip in Eq. (|27a[) is negligible such that 
we can effectively use the "raw recipe" ((8]). Laser pulse 
length, wave packet, and measuring point have been cho- 
sen such that the wave packet runs through the measur- 
ing point at a time where the laser is still fully active. 
It becomes obvious that the result from the wave func- 
tion in x-gauge is unphysical. The source of the prob- 
lem lies in the contribution exp (ip(t)zM)- As soon as 
the momentum p = p(t) moves even slightly, the pos- 
sibly large zm can amplify such a small oscillation and 
induce dramatic phase oscillations which, in turn, pro- 
duce a large contribution to the PES. This contribution, 
however, must be unphysical because it sensitively de- 
pends on the choice of the measuring point. Clearly, the 
v-gauge is the preferred choice for the evaluation of PES. 
This was already expected from the analytical consider- 
ations in section III C 31 We will henceforth exclusively 



use the u-gaugc. Of course, a numerical solution of the 
Schrodinger equation in coordinate space is often much 
simpler in x-gauge. In such a case, one still can use the 
x-gauge for the solution and then use the gauge transfor- 
mation (fl"7|) to bring this into u-gauge. This is the path 
actually followed in section IIII Bl 

Results for higher intensities are collected in figure [3] 
Here we use w-gauge throughout and compare the "raw" 
recipe (JSJ) with the PA recipe (|27a[) . The lowest panel 
shows the laser profile (|13b[) together with the squared 
wave function \ip(zM,t)\ 2 at the measuring point. This 
signals a critical situation where the laser is fully active 
while the wave packet is passing by the measuring point . 
The three upper panels show results for different inten- 
sities. The "raw" recipe still works acceptably well for 
the moderate intensity I = 10 12 W/cm 2 , but becomes 
grossly misleading for higher intensities. The generalized 
recipe (|27ap visibly improves the performance. The re- 
sults become reliable up to / = 10 13 W/cm 2 and remain 
somehow qualitatively correct for the highest intensity. 
It is thus much preferably to use the phase augmented 
form (|27a| for the evaluation of the PES. 

Figure [4] shows results from a configuration where 
the measuring point has been moved farther away to 
zm = 2000 which decouples the laser pulse from the wave 
packets signal. The case of highest intensity / = 10 14 
W/cm 2 is shown. As expected, the "raw" recipe yields 
the same result as the PA one and both agree nicely with 
the exact result. Moreover, the x-gauge (not shown here) 
yields precisely the same result as the w-gauge. Thus 
all distinctions and considerations are unnecessary in the 
case that laser pulse and particle flow do not overlap. 
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FIG. 3: Photoelectron spectra from wave packet analysis. 
Lower panel: The laser pulse envelope (dashed line) and the 
probability density \ip(zM,t)\ 2 at the measuring point (full 
line). Upper panels: PES analysis at zm m ^-gauge from 
"raw" (green or light gray curve) and with phase augmented 
(PA) signal (brown or dark gray line) for three laser intensi- 
ties as indicated. For comparison, the momentum distribution 
of the wave packet as given in Eq. (|37[) is also shown (black 
curve) . 
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FIG. 4: As in Figure [3] but for a situation where the laser 
pulse does not overlap with the wave packet at the measuring 
point (see bottom panel). Only the case of strongest laser 
intensity is shown. 



B. Realistic test cases 

As a first realistic test case, we come back to the in- 
troductory example of Figure [TJ namely the cluster Na^j" 
modeled with spherical jellium background and treated 
by time-dependent density functional theory using the 
energy functional of (25|. We take the laser with a typ- 
ical infrared frequency and two rather large intensities. 
The axial symmetry of this and of the following example 
allows us to perform the calculations in cylindrical coor- 
dinates. Figure [5] shows the results for the "raw" and the 
PA recipes. The "raw" PES in the upper panel repeats 
the case of figure [S] showing the obnoxious high energy 
shoulder. The PA PES makes a dramatic difference. It 
shows a reasonable, almost monotonous, decrease of the 
envelope. Only at about 2.5 Ry, the decrease turns into 
a rather weak slope which may be unrealistic as we come 
here into a region of very low yield where unwanted back- 
ground may spoil the analysis. The next lower intensity 
(lower panel) already shows reasonable pattern with the 
"raw" recipe. The PA recipe brings some improvement 
as it removes the glimpse of a shoulder at about 1.1 Ry. 
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FIG. 6: Ionization properties for the C4 chain under the 
influence of a laser pulse having frequency aji as = 0.1 Ry, 
pulse length T pu i sc = 12 fs, and intensity as indicated. The 
laser field was given effectively in w-gauge. The "raw" results 
(green of light gray curve) are given by Eqs. ((HJ, while the 
PA results (brown or dark gray line) are complemented by 
the phase factor (|27a|) . 



FIG. 5: Ionization properties for Nag" with jellium back- 
ground under the influence of a laser pulse having frequency 
'j-'ias = 0.1 Ry, pulse length T pu i sc = 12 fs, and intensity as 
indicated, computed in a cylindrical box of 176x88 Oq with 
spherical absorbing bounds covering at least 16 grid points. 
The laser field was given effectively in u-gauge. The "raw" re- 
sults (green or light gray curves) are given by Eqs. ((8}, while 
the PA results (brown or dark gray curves) are complemented 
by the phase factor (|27a|) . 



As already observed in the analytic case of a Gaussian 
wave packet, lower intensities perform well already with 
the "raw" recipe and the difference brought in from the 
PA recipe is negligible. 

Figure |5] shows a next test case results for the C4 chain, 
treated with non-local pseudo-potentials of Goedecker 
type [lH and, again, the electronic energy-density func- 
tional of [25| . The unphysical shoulders obtained with 
the "raw" recipe are even more developed than in case 
of Na^~ and their successful removal by the PA recipe 
is impressive. This is a clear demonstration of the gain 
achieved by the PA recipe (|27ap . 

IV. CONCLUSION 

We have investigated the evaluation of photo-electron 
spectra (PES) by sampling the time evolution of the wave 
function at a fixed measuring point. The schemes were 
first analyzed using solvable models. The traditional 
scheme was developed from the idea of a freely propa- 



gating outgoing wave. This requires negligible potentials 
at the measuring point which, in turn, sets a limit to 
an acceptable laser field strength. In order to extend 
the applicability of the scheme for evaluating PES to 
stronger fields, we have considered a model of a free par- 
ticle plus an external laser field in dipole approximation. 
This model is still analytically solvable and the analyt- 
ical solution allows us to deduce a more general scheme 
for the PES which consists in augmenting the collected 
wave function by an appropriate phase factor account- 
ing for the time dependent laser field. We have also 
investigated the effect of gauge transformations on the 
results and pondered the question of the most appropri- 
ate choice of gauge. We have found that the appropriate 
gauge for evaluating PES is clearly the v-g&uge (velocity 
gauge). Still, the most practical way to produce the wave 
functions in w-gauge is to compute them first in x-gaugc 
(space gauge) and to apply then the appropriate gauge 
transformation. 

The augmented scheme has been tested in an analyti- 
cal model of Gaussian wave packets and in a few realistic 
examples. A major result is that the "raw" recipe for 
PES is valid for low and moderate laser fields (about 
10 12 W/cm 2 for the test case Na^"). It remains valid in 
case that the particle flow arrives at the measuring point 
after the laser pulse has died out (in this case for all in- 
tensities). The phase augmented evaluation of PES was 
shown to considerably extend the range of applicability. 
The gain was particularly dramatic for the example of 
the C4 chain. With the augmented recipe for computing 
PES we have been able to check laser intensities up to 
10 14 W/cm 2 for our realistic test cases which all involved 
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a low laser frequency of 0.1 Ry. Higher laser frequencies 
reduce the effective intensity (Keldysh criterion). Thus 
even higher intensities may be used for higher frequen- 
cies. 

After the successful tests shown here, the augmented 
recipe for evaluating PES is ready for use in more 
demanding situations as they are molecular systems 
with large ionization potentials for which experimental 
data already exist as, e.g., N2, Cgo or typical organic 
molecules. 



Appendix A: Solution of the Schrodinger equation 
for the wave packet model 

We provide the analytical solution of the wave packet 
propagation in x-gaugc. This is simpler and the i>-gaugc 
is regained easily by the phase transformation (|17|) . 

The starting point is the Schrodinger equation as given 
in Eq. (|18a[) . The ansatz for the solution is given by 
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with fj,(t) = [Aq + i£t. The other time-dependent pa- 
rameters are determined substituting (|A1|) in the time- 
dependent Schrodinger equation and comparing term by 
term. In order to achieve this we first build the necessary 
derivatives: 
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Identifying term by term, we obtain the following equa- 
tions for the parameters 
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